Observation of Stiickelberg oscillations in accelerated optical lattices 
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We report the experimental observation of Stiickelberg oscillations of matter waves in optical 
lattices. Extending previous work on Landau-Zener tunneling of Bose-Einstein condensates in optical 
lattices, we study the effects of the accumulated phase between two successive crossings of the 
Brillouin zone edge. Our results agree well with a simple model for multiple Landau-Zener tunneling 
events taking into account the band structure of the optical lattice. 
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Stiickelberg oscillations arise when a quantum system 
is forced more than once through an avoided crossing 
of two of its energy levels [l-Sij. Each passage through 
the avoided crossing leads to Landau-Zener tunneling in 
which the initial wavefunction is coherently split between 
the two energy levels, and in the next passage the phase 
accumulated between the energy levels leads to quantum 
interference. As a result, the populations in the two en- 
ergy levels after several Landau-Zener passages are very 
sensitive to the dynamics between the passages. While 
Stiickelberg oscillations were first predicted more than 
eighty years ago for atomic collisions, they have received 
renewed interest in recent years as a tool for the spec- 
troscopic investigation [1] of, e.g., superconducting arti- 
ficial atoms molecule formation in ultra-cold gases 
[1], and Rydberg excitations [S]. Here we report results 
on the realization of Stiickelberg oscillations using Bose- 
Einstein condensates in accelerated optical lattices. From 
the measurements of the frequency and damping time of 
the oscillations we are able to deduce the distance be- 
tween the energy bands of the optical lattice and the 
momentum spread of the condensate. 

In our experiments we investigate the time evolution 
of Bose condensates loaded into a spatially periodic po- 
tential subjected to an additional (inertial) static force, 
as described in detail in [l^. The dynamics of ultracold 
atoms in a tilted optical lattice can be described by the 
well-known Wannier-Stark Hamiltonian [ll| determining 
the atomic evolution in an optical lattice of depth Vq and 
spatial period c^l = Al/2 determined by the wavelength 
Al of the laser creating the periodic potential. The char- 
acteristic energy scale of the system is the recoil energy 
Ei-cc — tt^Tt? /2Md^, where M is the atomic mass. The 
band-structure representation of the energy En {q) of the 
n-th band versus the quasimomentum q is represented 
in Fig. [TJ The observation of Stiickelberg oscillations 
is based on a sweep of the quasimomentum q given by 
q(t) = Ft induced by a force F (measured here in units 
of E^cc/dL) acting on the atoms in the rest frame of the 
lattice. 

Our experiments can be modeled (neglecting the finite 
duration of the Landau-Zener tunneling event [lol, IHI) 
using a two-state ansatz for the two lowest energy bands 



EQ,i{q) in the optical lattice, described by the state vec- 



tors I Q I and 



The different steps of the ex- 



perimental protocol are then represented by matrices 

a 



acting on the general state vector |^ ^ j [Ic 
• LZ tunneling event: 



1lz(/^lz) = 



^/T^rp^e,-i{'pst-^n) 



^/Y-775^e*(0st-^/2 



where Plz — exp(— 7r/7) is the Landau-Zener (LZ) 
tunneling probability with 7 = 'i2F di^E^^c / {t^Vq) , 
and 0st = 7r/4 - [I + ln(2^\[(2-i)-^ + argT[l - 
(27)^1^] the Stokes phase 

• acceleration from initial quasimomentum qini to fi- 
nal quasimomentum ggn: 



Macc(gini,gfin) 



-j*o(9ini,gtin) 







,-i$i(gi„i,gfi„) ) ' 



where, using the linear dependence of q on time, 

*n(*ni, gfin) = ^ En{q)dq. 

• waiting for time iwait at quasimomentum q: 



'wait {q, t-wait) 



g-iBo(g)t„ait/'i Q 

Q g-i-Ei(9)twait/?i 



A general experiment of crossing the zone edge at 
q — 0.5 qB, accelerating to qtin, waiting for time twait, 
accelerating back to q — 0.5 qs and crossing the zone 
edge a second time (see Fig. 1) can then be described by 
the total transfer matrix 



wait '^'acc Mlz, 



(1) 



where the arguments of the matrices have been omitted 
for brevity. If the system is initially in the lowest energy 
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FIG. 1: Experimental protocol for observing Stiickelberg os- 
cillations in accelerated optical lattices. 



FIG. 2: Stiickelberg oscillations after two crossings of the 
zone edge in an optical lattice. Shown here is the survival 
probability in the lowest energy band Po for Vb ~ lAErcc and 
F = 1.197 Ercc/dh as a function of qa-n; for these experimental 
parameters the Stokes phase is <j!>st ~ — vr/3. The dashed line 
is a numerical simulation assuming a momentum spread of the 
condensate of Ag = COSgB- The data points are the results 
of single experiments and hence no error bar is reported; shot- 
to-shot fluctuations for constant experimental parameters are 
around 10%. 



band of the lattice, the probability Pq of finding it in that 
state after the above protocol is 

Po = 4Plz(Plz - 1) sin2 + cj)^^ + 1, (2) 

where = [($o - + (-^o - i?i)iwait/?j] denotes the 
total relative phase between the two bands accumulated 
during the acceleration and waiting parts of the protocol. 

Our experimental setup is described in detail in [lo| . 
Briefly, we create Rb condensates containing up to 10^ 
atoms in a crossed dipole trap and subsequently load 
them into a one-dimensional optical lattice formed by two 
linearly polarized counter-propagating gaussian beams 
far-detuned from the Rb atomic resonance. The dipole 
trap beam perpendicular to the lattice direction is then 
switched off and the lattice is accelerated by applying 
appropriate frequency chirps on the radio-frequency gen- 
erators powering the two acousto-optic modulators used 
to control the power in the lattice beams. Finally, the 
dipole trap beam and the lattice beams are switched off 
and an absorption image is taken after 23 ms of time-of- 
flight, which allows us to deduce the survival probability 
in the ground state Pq = /^tot from the relative num- 
ber of atoms in the p = momentum class. 

Stiickelberg oscillations are then observed by first load- 
ing the condensate adiabatically into the lowest energy 



band of a lattice (Vq = 1.4Proc) with zero quasimomen- 
tum (stationary lattice), accelerating the lattice across 
the edge of the Brillouin zone at q = 0.5 (where 
(7b — lirfi/di^ is the Bloch momentum), reversing the sign 
of the acceleration, crossing the edge of the Brillouin zone 
a second time and finally stopping at zero quasimomen- 
tum. The results of such an experiment with iwait = 
(i.e. Mwait = I) and varying qf^^^ are shown in Fig. 2. 
Since the energy separation between the bands varies (ap- 
proximately) linearly with g^n, the relative phase is 
a quadratic function of qf^r^ and hence the frequency of 
the oscillations varies, leading to 'chirped' Stiickelberg 
oscillations. Beyond quasimomentum q = in Fig. 1 
our simple two-level model is, strictly speaking, no longer 
valid as avoided crossings with higher bands will lead to 
a loss of population from the excited band. Since the 
tunneling probability to those higher bands (and back) 
is close to unity, however, the dominant correction to the 
two-level approximation will be in the phase factor. 

In spite of the simplifying assumptions, the agreement 
with our simple model (using Mathieu-function solutions 
for the band structure of the lattice) is still very good. 
From the total acceleration time tacc = 0.4 ms between 
successive crossings of the zone edge for the maximum 
gfin, the damping time of the Stiickelberg oscillations of 
Fig. 2 is found to be around t = 0.35 ms. We can account 
reasonably well for this damping by assuming a momen- 
tum spread of the condensate of Ag = 0.03gB in our sim- 



3 



ulations, which in practice is partly due to the intrinsic 
momentum width of the condensate and partly to de- 
phasing mechanisms such as dynamical instabilities that 
occur during the acceleration phase The momen- 

tum spread leads to an accumulation of relative phase 
of different parts of the condensate at a different rate 
depending on their position in the Brillouin zone. We 
can, therefore, connect the momentum spread to the 
spread in the accumulated relative phase A0 in the fol- 
lowing way: Assuming that the two lowest energy bands 
for a shallow lattice correspond roughly to the free par- 
ticle dispersion relations En{q) = {q — nqB)'^ / {2M) (ex- 
cept around q = 0.5 qs), we find dA(p/dq = 8-Erec7tot/?i 
with Ttot — 2tacc + twait • For a given momentum spread 
Aq, the damping time t of the Stiickelberg oscillations 
is then r = h/ {^E.^dAq/ qs)), giving Aq/qB = 0.02 for 
T = 0.35 ms, in good agreement with the simulations. 
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FIG. 3: Survival probability in the lowest energy band Po 
for Vo = 1.1 Ercc and F — 1.197 iJrcc/dL as a function of the 
waiting time at q^n = 0.719 gs. 

We now fix qn^ = 0.719 qs and modify the above proto- 
col by inserting a variable waiting time between the two 
acceleration stages. As can be seen in Fig. 3, this gives 
rise to Stiickelberg oscillations with constant frequency 
since A(f> now varies linearly with iwait- The measured 
frequency of the oscillations of around 6.7 kHz is in good 
agreement with the expected band separation of 7.07 kHz 
at g = 0.719 qs- As in Fig. 2, the oscillations are damped 
on a timescale of ~ 0.4 ms. 

Finally, we perform an experiment in which the lattice 
depth Vo is varied keeping qn^ — 0.757 qs and t^ait ~ 
constant. In this case, both Mlz (through Plz) and Mace 
(through the band structure and hence Eo_i{q)) are var- 
ied and the resulting variation in Pq depends sensitively 
on the exact parameters of the system. In the two ex- 
treme cases of very small and very large lattice depth Vq , 
however, the experiment has an intuitively predictable 



FIG. 4: Survival probability in the lowest energy band Po as 
a function of Vo for F = 1.197 Srcc/dz,, a fixed final quasi- 
momentum gfln = Q.Tal qs and fwait = 0. The Stokes phase 
varies between — 7r/2 for Vo ~ 0.5 iSrcc and 7r/9 for Vb « 4 iSrcc- 



Eq. 2). For small Vb, the Landau-Zener tunneling prob- 
ability Plz ~ 1 and hence complete tunneling will occur 
in both Landau-Zener passages, so the system ends up 
in its initial state after the above protocol. For large Vq^ 
we also expect the system be in the ground state energy 
band after two Landau-Zener passages, but now this is 
because P^z ~ and hence tunneling into the other band 
is strongly suppressed, so the system stays adiabatically 
in the ground state at all times. 

Fig. 4 shows the results of the above experiment. As 
expected, for very large and very small values of Vb, is 
close to 1, whereas in between Pq reaches a minimum of 
« 0.3 for Vb ~ 1.4 i^rcc- We note here that, as can be seen 
from Eq. 2, the exact behaviour of Pq depends crucially 
on how (which includes the Stokes phase) depends 
on Vb- Indeed, if no relative phase were accumulated 
between the two crossings of the zone edge, we would 
expect Po = 1 independently of the value of Plz- 

In summary, we have observed Stiickelberg oscillations 
of Bose-Einstein condensates in optical lattices using dif- 
ferent experimental protocols. The sensitivity of those 
oscillations to accumulated phases in the energy bands 
suggests the use of Stiickelberg protocols for studying 
dephasing and related phenomena in optical lattices, for 
instance in a superlattice [11]. Our methods for the de- 
termination of dephasing and energy splitting may be 
also be applied to other systems that can be described 
by a two-state model. 
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